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Abstract. Possible scenarios based on available experimental data and phenomenological 
knowledge of the GeV energy region are extended to the TeV energy region in the framework 
of the weighted superposition mechanism of soft and semi-hard events. KNO scaling violations, 
forward-backward multiplicity correlations, H q vs. q oscillations and shoulder structures are 
discussed. 



1. Introduction 

Given the current difficulty of performing QCD calculations in the realm of multiparticle 
dynamics, one possible alternative path consists in studying the features of hadronic final 
states as obtained in collision experiments. The main motivation is the conviction that the 
complex structures which we observe might very well be simple at the origin, and that such 
initial simplicity manifests itself in terms of regularities in final particle multiplicity distributions 
(MD's). 

Indeed, the Pascal regularity ^ appeared very soon: that final charged particle MD's were 
well described by the negative binomial (NB), also known as Pascal, distribution was discovered 
first in cosmic rays observations [2j in the '60, and later confirmed in the 5-100 GeV energy 
range in pp collisions and up to 40 GeV in e + e~ annihilations, as well as in other types of 
collisions (see, e.g., 

The NB (Pascal) distribution is a two-parameter distribution: 

^ B -'(«.*) = * ( * +1) "J* +n "" (a ^ m 

where n is the average multiplicity and 1/k measures the deviation of the variance D 2 = (n 2 ) —n 2 
from the Poisson shape: 

l/k + 1/n = D 2 /n 2 . (2) 

In fact, for the Poisson distribution D 2 = n, i.e., k — > oo, and for the geometric distribution 
D 2 = n + n 2 , i.e., k = 1. 

That different reactions showed the same, approximate, regularity was worth of interpretation 
attempts. Probably the most successful of such attempt has been clan structure analysis (I]. 
'Clan' (recalling the Scottish sense of the word) refers to a group of particles of common ancestry: 
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Figure 1. The shoulder structure visible in this charged particle MD at 900 GeV cm. energy 
from the UA5 Collaboration ,5i is well described by the weighted superposition of two NB 
(Pascal) MD (dashed and dotted lines). 

clans are by definition independently produced in a number which follows the Poisson MD, each 
clan contains at least one particle by assumption and all correlations are exhausted within 
each clan. Clan ancestors, after their production, generate additional particles via cascading, 
according to a logarithmic multiplicity distribution. 

The two relevant parameters in clan structure analysis are the average number of clans, N, 
and the average number of particles per clan, n c . They are linked to the standard NB (Pascal) 
parameters n and k by the following non-trivial relations 

N = fcln(l + n/fc) and n c = n/N. (3) 

Clan structure analysis reveals new interesting properties when applied to above mentioned 
collisions. In particular it is shown that in e + e~ there are more clans than pp collisions, whereas 
each clan is much smaller. In addition, clans in central rapidity intervals are larger than in more 
peripheral intervals. The deep inelastic case is intermediate between the previous ones: clans 
are less numerous than in e + e~ annihilation and hadronic in character, but the clan size tends 
to be leptonic. 

As soon as the cm. energy increased, more detailed structures appeared in the data: the MD 
was seen to have a 'shoulder structure' in the intermediate multiplicity range, first in pp collisions 
|H] and later also in e + e~ annihilation J . In both cases an explanation was found in terms of 
the weighted superposition of different classes of events: events with and without mini-jets in the 
case of pp collisions (called 'semi-hard' and 'soft', respectively; see Figure^), events with a fixed 
number of jets in e + e~ annihilation |S]. It should be stressed here that the MD of each class 
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Figure 2. Results of the weighted superposition model for ne(rtf) vs. uf compared to 
experimental data in full phase-space at 63 GeV (a) and in the pseudo-rapidity interval |r/| < 4 
at 900 GeV (b) HHJ. 



was again well described by the Pascal distribution (which could not, of course, by itself alone 
reproduce the shoulder structure). The degree of precision of such a description is witnessed 
by the fact that not only the overall shape is visibly reproduced, but also the quasi-oscillatory 
behaviour of the so-called H q moments (defined as the ratio of factorial cumulant moments of 
order q to factorial moments of the same order, H q = K q /F q ), which show sign changes when 
plotted ciS ct function of the order q (§| . 

Finally, a third observable was brought into play, with different characteristic behaviour in 
e + e~ annihilation and pp collisions, but which can be understood in terms of the Pascal regularity 
and clan structure analysis: the strength of forward-backward multiplicity correlations jlOj . In 
particular, it was shown that in e + e~ annihilation the superposition mechanism is sufficient by 
itself to account for the small amount of forward-backward multiplicity correlations (FBMC) 
present in the data, consistently with the idea of a large number of small clans; on the other 
hand, the superposition was not enough in pp collisions (Figure EJ): it is required that particle 
be produced in sizeable groups (clans!) which are not localised in phase-space, i.e., clans must 
produce particles in both hemispheres (in other words, there must be a non-zero 'leakage': there 
are particles going backwards from forward-emitted clans, and vice versa.) 

2. Extrapolations to higher energies 

We will now abandon the parallelism with e + e~ annihilation, and in the following concentrate 
on pp and pp collisions. The knowledge from the GeV energy range, summarised briefly in 
the previous section, will now be employed to extrapolate the MD's in the TeV energy region 
j 1 11 112j . i.e., the region covered by Tevatron and LHC colliders. 
The main equation of the extrapolation is the following: 

^ total) = «sof tJ Pi Pascal) (ni, h) + (1 - a soft )Pi Pascal )(n 2 , k 2 ), (4) 

where a S oft is the fraction of 'soft' events in the total sample. An event is declared 'soft' if it 
contains no mini-jets: there is some ambiguity as to the threshold which defines a mini-jet, e.g., 
a calorimetric tower of a few GeV, but a small variation of this threshold does not result in 
a wild variation of each class' parameters. On the other hand, a proper definition of what is 
"hard" and what is "soft" is an extremely interesting task on itself, but it goes well beyond the 
scope of this paper. 
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Figure 3. Extrapolation of the average 
multiplicity. 
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Figure 4. Extrapolation of the weight 
parameter. 



The basic assumptions to be used in the extrapolation are as follows (see Figure 0J): 

(i) the average charged multiplicity in the total sample grows with the logarithm of the cm. 
energy (best fit to available data); 

(ii) in agreement with the findings at ISR, where one component gives an adequate description 
of the data, the growth of the average charged multiplicity in the soft sample is proportional 
to In s. 

(hi) in agreement with UAl findings at SpS, the average multiplicity in the semi-hard sample 
is twice as large as in the soft one, but a very small term quadratic in In s could be added 
as a correction (the two possibilities are called A and B in the figures). Accordingly, most, 
if not all, of the In 2 s behaviour of fttotal comes from the increase in the mini-jet production 
cross section. 

One now has enough information to compute the expected energy dependence of the weight 

a so ft(s) = 2 - n tota i(s)/n soft (s), (5) 

which is shown in Figure |1] 

We come now to the behaviour of k. k so f t was found to be constant in the GeV region by 
the UA5 collaboration; together with the mentioned behaviour of n so ft this implies that KNO 
scaling behaviour is valid for the soft component in the TeV region: 

^soft/^soft ~ 0- 14 « constant. (6) 

We stay with this assumption on fc so ft. 

As the data on the semi-hard component are scarce, it was decided to explore three different 
scenarios for the energy behaviour of & S emi-hard> summarised, together with the average number 
of clans and the average number of particles per clan, in Figures 15171 

Scenario 1 (Figure^) 

In this scenario we assume that KNO scaling holds also for the semi-hard component, i.e., that 
fcsemi-hard is constant. This leads to a slow increase of the average number of clans and of the 
average number of particles per clan. 
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Figure 5. Interpolated and extrapolated NB (Pascal) parameters in scenario 1. 



Scenario 2 (Figure^ 

Here we test the assumption of maximum KNO scaling violation, by assuming that the growth 
fctotai — 0.0512 In — 0.08 continues in the TeV region. One notices that this scenario implies a 
fast decrease with increasing energy of the average number of clans, a fact which is quite unusual 
and interesting, and is explored further in these proceedings |13j . 

Scenario 3 (Figure^ 

In this scenario we transport the QCD-predicted (at leading log level) behaviour to fit the low 

energy pp data: k~ e l ai _ haid = 0.38— \J 0.42/^(^/10). It corresponds to an increase, slower than 
in scenario 2, toward an asymptotic constant value, i.e., toward a MD which is asymptotically 
KNO scaling. From the point of view of clan analysis, this scenario is intermediate between 1 
and 2, and shows again a decrease in N with increasing y/s (N will start to increase again in 
the KNO scaling regime, due to the increase of n.) 



3. Extension to small rapidity intervals 

In going from full phase-space (FPS) to pseudo-rapidity (??) intervals, our main concern is 
to be consistent with the scenarios explored in FPS and extend them. It should be pointed 
out that by assuming only a longitudinal growth of phase space and constant height of the 
rapidity plateau with cm. energy for semi-hard events, as done in Ref. J2j, CDF data Jl] in 
pseudo-rapidity intervals are underestimated. These data are well described by allowing a In 2 s 
growth of the total rapidity plateau: from this consideration one deduces a more appropriate 
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Figure 6. Interpolated and extrapolated NB (Pascal) parameters in scenario 2. 



growth of the semi-hard plateau height; the constraint is that ra S emi-hard i n full phase space 
follows a logarithmic growth with sfs as discussed above (curves B). Predicted charged particle 
multiplicity distributions in the three scenarios of the two-component model for the interval 
\r}\ < 1 at Tevatron and LHC energies, calculated for the last mentioned case, are shown in 
Fig. IH1 (notice that a direct comparison with CDF data, in view of the relatively large value of 
their resolution p T > 0.4 GeV/c, is questionable). If this behaviour for the semi-hard component 
will be confirmed by data one should conclude that semi-hard events populate mainly the central 
rapidity region giving an important contribution to the increase of charged particle density in 
central rapidity intervals. 

4. Fermilab results 

Two sets of data are available from Tevatron: the E735 Collaboration ^5] gives full phase-space 
results, from data measured in \r]\ < 3.25 and p T > 0.2 GeV/c and then extended via a Monte 
Carlo program, which do not completely agree with those obtained at comparable energies at 
the SpS collider, as shown in Figure El (left panel); Tevatron data are more precise than SpS data 
at larger multiplicities (they have larger statistics and extend to larger multiplicities than UA5 
data), but much less precise at low multiplicity. Both sets of data show a shoulder structure, 
but the Tevatron MD is somewhat wider. As seen from the same Figure El (right panel), scenario 
2 is the one that come closer to these data. 

The second set of data comes from CDF ^Hl for charged particles in \rj\ < 1 and p T > 
0.4 GeV/c: it was found that by subdividing the minimum bias sample into two groups, 
characterised respectively by the absence ('soft' events) or the presence ('hard' events) of mini- 
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Figure 7. Interpolated and extrapolated NB (Pascal) parameters in scenario 3. 




Figure 8. E735 results do not quite agree with UA5 ones (left panel; data from the two 
experiments which were taken at nearly the same energy are here rescaled by the same factor), 
but are close to scenario 2's predictions (right panel). 
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Figure 9. Predictions for the multiplicity distributions in \rj\ < 1 at 1800 and 14000 GeV in 
our scenarios for pp collisions. 



Table 1. NB (Pascal) and clan structure 
analysis parameters for pp collisions at 14 TeV 
in the two-component model. 
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Table 2. Forward-backward multiplicity cor- 
relations strength parameter for pp collisions 
at 14 TeV in the two-component model. 



FB corr. strength 
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Figure 10. CDF results on MD's at 1800 GeV: the soft component satisfies KNO scaling, the 
hard one does not [To] . 



jets, interesting features of the reaction can be investigated. More precisely, a 'hard' event has 
been defined as an event with at least one calorimeter cluster in \rj\ < 2.4, a cluster being defined 
as a seed calorimeter tower with at least 1 GeV transverse energy Et plus at least one contiguous 
tower with Et > 0.1 GeV. A subdivision which is interesting per se and can be tested at 14 TeV. 
In summary, the soft component is found to satisfy KNO scaling (as expected in our scenarios), 
while the hard one does not (in disagreement with scenario 1, but in agreement with scenarios 
2 and 3); also the {p ) distribution scales at fixed multiplicity in the soft component and not 
in the hard one; the dispersion of {p ) vs. the inverse of the multiplicity is compatible with an 
extrapolation to as n — > 00 in the soft component but not in the hard one, indicating in this 
limit a lack of correlations in the soft component. 

5. Conclusions 

The weighted superposition mechanism of two NB (Pascal) MD's describes not only MD and 
H q oscillations, but also forward-backward multiplicity correlations in pp and pp collisions and 
e + e~ annihilations. In hadronic collisions, the two components were found to correspond to soft 



events (without mini-jets) and to semi-hard events (with mini-jets) respectively. Based on this 
mechanism, the knowledge of the features of MD's up to 900 GeV cm. energy has been used 
to predict the characteristic behaviour expected in the TeV energy range: the soft component 
satisfies KNO scaling, while the semi-hard one violates it strongly. These predictions do not 
disagree with the successive findings at Tevatron. In terms of clan structure analysis, a peculiar 
behaviour was found: particles in the semi-hard component tend to aggregate forming few, 
large clans. This behaviour, and its consequences for the LHC experiments, will be investigated 
further in A. Giovannini's contribution to these proceedings |13j . 
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